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Abstract 

Motivated by the interest in searches for violation of CPT invariance, we study its possible effects 
in the flavour ratios of high-energy neutrinos coming from cosmic accelerators. In particular, we 
focus on the effect of an energy-independent new physics contribution to the neutrino flavour 
oscillation phase and explore whether it is observable in future detectors. Such a contribution 
could be related not only to CPT violation but also to a nonuniversal coupling of neutrinos to 
a torsion field. We conclude that this extra phase contribution only becomes observable, in the 
best case, at energies greater than 10 16 ' 5 GeV, which is about five orders of magnitude higher than 
the most energetic cosmological neutrinos to be detected in the near future. Therefore, if these 
effects are present only in the oscillation phase, they are going to be unobservable, unless a new 
mechanism or source capable to produce neutrinos of such energy were detected. 
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I. INTRODUCTION 



Experimental evidence has confirmed that flavour transitions are the solution to the for- 
mer so-called solar and atmospheric neutrino deficit problems [2], [J. Further evidence was 
provided by experiments performed with neutrinos generated in particle accelerators and 
nuclear reactors, such as KamLAND^] and K2K|5]. The mechanism responsible for these 
transitions requires neutrinos to be massive: the probability of a flavour transition is oscil- 
latory, with oscillation length X std = AirE / Am 2 , where E is the neutrino energy and Am 2 
is the difference of the squared masses of the different neutrino mass eigenstates. However, 
even though this mass-driven mechanism is the dominant one in the energy regimes that 
have been explored experimentally (MeV-TeV), there is still the possibility that alternative 
mechanisms contribute to the flavour transitions in a subdominant manner, which perhaps 
can manifest at higher energies. 

Although these alternative mechanisms involving new physics (NP) are able to produce 
flavour transitions, it is known that none of them can explain the combined data from 
atmospheric, solar, accelerator, and reactor neutrino experiments performed in the MeV-TeV 
range, unlike the pure Am 2 oscillation mechanismlg, M, M, M . Some of these alternatives! 101 

n □ 

are the violation of the equivalence principle (VEP), of Lorentz invar iance[ll|, [12j (VLI), 
of CPT invariance (VCPT), the non- universal coupling of neutrinos to a space-time torsion 
field (NUCQ), decoherence during the neutrino's trip, and non-standard interactions (NSI). 

Typically, these mechanisms result in oscillation lengths \ NP that have a different depen- 
dence on E, usually expressed as a power-law, A ~ E n , with the value of n depending on 
which mechanism is being considered: for instance, n = for VCPT and NUCQ and n = 1 
for VEP and VLI, while with n = — 1 the standard Am 2 oscillations are recovered. Atmo- 
spheric events from Super-Kamiokande (SK)^] were used to find the value n = —0.9 ± 0.4 
at 90% C.L., thus confirming the dominance of the Am 2 oscillation mechanism and forcing 
any other mechanisms to be subdominant, at least within the energy range and pathlength 
considered in said analysis. 

So far, searches for NP effects in neutrino oscillations have been limited to energies 
ranging from a few MeV to a few Gevjfj, Q, 9] and have turned out negative. How- 
ever, proposals for analyses of atmospheric neutrinos with energies of up to 10 4 TeV in 
second-generation neutrino detectors such as IceCubefl3| and ANTARES[14j are being con- 
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sidered. Due to the energy dependence of the oscillation lengths, the oscillation phases scale 
as (27rL/A 7VP ) / (27r/A s * d ) ~ E 1 ^ 11 , that is, the relative dominance of the NP contribution 
grows with the neutrino energy provided that n < 0, so that the observation of very energetic 
neutrinos -such as the ones expected from presumed cosmic accelerators like active galaxies 
and gamma ray bursts- would offer a means to establish whether the Am 2 oscillation mech- 
anism is still the dominant one at high energies or to otherwise set stronger bounds on the 
NP parameters. 

In this work, we have introduced the aforementioned new physics through the use of a 
modified dispersion relation, and focused our analysis on the case of n — 0, corresponding, 
as we will see, to an energy-independent NP contribution to the neutrino flavour oscillation 
phase. We have calculated the proportion of each flavour arriving at Earth from distant 
cosmic accelerators and explored how it is affected by the parameters that control the new 
physics, and whether these effects are observable at all. 

The outline of the paper is as follows. Section [Til describes how the NP arising from 
a modified dispersion relation affects the flavour-transition probability for neutrinos that 
travel cosmological distances. In Section II I II we explore the case of an energy-independent 
NP contribution and its effects on the flavour ratios. Finally, in Section HV] we present our 
conclusions. 



II. FLAVOUR-TRANSITION PROBABILITY IN THE PRESENCE OF A MOD- 
IFIED DISPERSION RELATION 

n 

The NP effects can modify flavour transitions in two ways|14l]: by transforming both the 
oscillation length and the neutrino mixing angles or by altering only the first. The former 
case occurs, for instance, when considering the low energy phenomenological model of string 
theory, known as 'Standard Model Extension' [lo| and has been examined using SK and K2K 
datap, Q. The second case can be achieved by considering a modified dispersion relation 
which departs from the well-known formula E 2 = p 2 + m 2 . Because we wish to explore 
whether solely effects on the phase are observable at high energies, we follow this second 



alternative and consider the following modified dispersion relation 15], which allows us to 
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study the contributions of NP effects in a model-independent way: 

(E \ a 
] = p 2 + m 2 + rip 2 E a , (1) 
mpj 

where mp ~ 10 19 GeV is the Planck mass, = rjmp is an adimensional parameter that 
controls the strength of the NP effects and, following the literature, a has been chosen to be 
of integer value. Such a dispersion relation assumes that the scale of NP effects is the Planck 
scale where, according to theories of quantum gravity, space-time might become "foamy". 
Eq. (TjQ) was recently[lJ] used to predict the sensitivity of the ANTARES neutrino telescope 
to NP effects in the high-energy atmospheric neutrino flux. 

We now derive the flavour-transition probability in the presence of NP effects, for neutri- 
nos that propagate over a cosmological distance. Flavour transitions arise as a consequence 
of the fact that flavour eigenstates |z/ a ) (a = e, /z, r) are not also mass eigenstates |z/j) 
(i = 1,2,3), but rather a linear combination of them, i.e. \v a ) = Y^=i ^mN)j w hh 
elements of the neutrino mixing matrix. 

Using the standard dispersion relation, it is a common procedure to derive an approximate 
expression for the momentum of the z-th neutrino mass eigenstate, 

2 

mf 



p l = ^E 2 -m 2 ^E-^, (2) 

where rrii is the mass of the neutrino and E is its energy, and are such that, at the energies 
that we have considered, m, <C E. From this equation we obtain the usual expression for 
the momenta difference: 

Am?- 

A Pij = p • - Pi = — — . (3) 

n 

In accordance with the latest bounds obtained from global fits[16j, we have set the three 
mixing angles that parametrise U to sin ($12 ) = VolOA, 0i3 = and 9 23 = tt/4. The mass- 
squared differences have been set to Am^ = 8.0 x 10~ 5 eV 2 and Am^ = 2.5 x 10 -3 eV 2 , and 
we have assumed a normal mass hierarchy (i.e. m 3 > mi), so that Am%i = Am^ + Ara2i. 
The probability that a neutrino created with flavour a is detected as having flavour (3 after 
having propagated a distance L in vacuum is given byjl?)] 



(4) 



Pu^u, (E, L)=5 aP -A-£ Re ( jf) sin 2 , 

i>j ^ ' 

where jff = U^JJ^U^U^. Since 13 = 0, U is a real matrix, independent of the CP 
violation phase, 5cp- 



In the case of the modified dispersion relation of Eq. ([TJ) we can also find an expression 
for the momenta difference. To first order in r)i, and discarding terms higher than second 
power in rrii or involving ^mf , we obtain 

mf T]iE n 

P ^ E -2E-—> (5) 

with n = a + 1, and hence 

Am 2 . An (n) E n 

where A^™ = r]^ — T]^ n \ for the NP mechanism with an E n energy dependence. Note 
that it is necessary that the rji have different values for different mass eigenstates in order to 
have a nonzero NP contribution to the momenta difference. The corresponding oscillation 
probability is Eq. (j4]) with Apy — > Ap^-; hence, the NP affects solely the oscillation phase, 
but not its amplitude. 

Since L>1 for high-energy astrophysical neutrinos, sin 2 (ApyL/2) is a rapidly oscillating 
function and so, due to the limited energy resolution of neutrino telescopes, the average 
flavour-transition probability is sometimes used instead, i.e. 

{Pu^ fi ) = Y t \U ai \ 2 \U P i\ 2 . (7) 

i 

Let (P Va -+ V0 ) std be the standard average probability, that is, when there are no NP effects 
present. If the extra term in ApY,- has n > 0, its effect will be that, at high energies, 
P Va _> v will oscillate even more rapidly with energy, but still around the same mean value 
(Pu a -wp) ■ If n < 0, the oscillations will continue up to high energies, also around the 
same mean, and at a high enough value, when Ap.^ = Amf-/2E — > 0, the probability will 
tend to zero. However, if n = 0, then the extra term in Apij is energy-independent and so 
when Ap,^ — > 0, P Va ^ Vfj becomes constant, but different from zero due to the existence of 
the extra term. Furthermore, when this happens, and depending on the values of the Ar^, 
it is in principle possible for the constant probability to be different from (P Ua _, U0 ) std . Such 
a nonzero, constant probability at high energies could therefore be interpreted as being due 
to the contribution from energy-independent new physics. We will focus on this possibility. 
Although our analysis of NP effects using Eq. (TjQ) is model- independent, A?/ ) takes a 



different form depending on the particular mechanism being considered 0, 



13| . In the energy- 



independent oscillation mechanism that we are focusing on, the extra contribution could be 
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due to VCPT through Lorentz invariance violation, in which case 

A^ 0) = h - h = 6 y , (8) 

i — I 

with hi the eigenvalues of the Lorentz- violating CPT-odd operator pi vi^^^y^. Alter- 
natively, the contribution could be due to NUCQ, and in this case|l8j] we would consider 
different couplings, fcj 7^ kj (for mass eigenstates i and j), to a torsion field Q, so that 

A vlf = Q{ki-kj) =Qhj. (9) 

n 

Strict bounds [6(1 have been set on the parameters that control the energy-independent NP 
mechanism using data from atmospheric and solar neutrinos, as well as SK and K2K, with 
energies up to about 1 TeV: 

621 < 1.6 x 1(T 21 GeV , 632 < 5.0 x 1(T 23 GeV . (10) 

Because the relative dominance of the NP energy-independent phase over the standard 
oscillation phase increases with neutrino energy, i.e. (Apy — Apy) /Ap^ ~ E, we would 
like to look at the most energetic neutrinos available. Hence, we will consider neutrinos 
originating at cosmic accelerators, such as active galactic nuclei, where it is presumed that 
they are created with energies of up to 10 11 GeV. Because the typical distance to these 
accelerators is in the order of hundreds of Mpc, we must include in the flavour-transition 
probability the effect of cosmological expansion. Hence, instead of the argument that appears 
in the sine of Eq. (J4j) , we define an (pij as follows: 



2E Q \t 10 E 



. v 5/3 /, \ 5/3' 



(U) 



where t{ and tf are the times at which the neutrino was produced and detected, respectively; 
t Q = 13.7 Gyr is the age of the Universe [20]]; and we have used the relation between the 
energy at detection (E ) and production epochs (E), in an adiabatically expanding universe, 
E (r) = E Q (t D /r) 2 ^ 3 = E (1 + z). Considering the detection time tf in the present epoch, 
tf = t , we obtain the accumulated phase 



Am?. [eV 2 ] 



^ v -„,~, ^[Gev] 



,- (E ,z) = 1.97 x 10 23 ""7_ L " n J l-(l + ^)" 5/2 , (12) 



where we have made use of the relation ti/t = (1 + z) 



-3/2 



By replacing the momenta difference Ap^ with Apij, we obtain, correspondingly, 



(frij (E a , z) = (fry (E , z) + 



1 - (1 + z 



n-3/2 



(E ,z)+$>{E ,z) , (13) 



with &j the contribution to the phase due to the NP effects. For n — 0, 



4° } (E , z) = 3.28 x 10 4 %- [GeV] 



41 j 



1 - (l + Z 



-3/2 



Hence, instead of the traditional expression in Eq. (TJJ for P Ua ^ U/3 , we will employ 

(S OI *) = <W - 4 £ Re (jf) sin 2 , 



(14) 



(15) 



where the explicit expression for 0^- = 0^ (E , z) is either of Eqs. (JT2"j) or ffT3]) . depending 
on which dispersion relation is being considered 251 ] . 



III. OBSERVABILITY OF THE NP EFFECTS IN THE HIGH-ENERGY NEU- 
TRINO FLAVOUR RATIOS 



Using the flavour-transition probability obtained in the previous section, Eq. (|15[) . we 
can calculate the ratio of neutrinos of each flavour to the total number of neutrinos that 
arrive at the detector from a source with redshift z. For a-flavoured neutrinos (a = e,fi,r) 
with energy E Q , this is 



(3=e,ji,T 



S 

1>P J 



(16) 



where is the ratio at the detector and is the ratio at the source. The latter is 
estimated assuming that neutrinos are secondaries of high-energy proton-proton or proton- 
photon collisions, which produce pions that decay into neutrinos and muons which decay into 



neutrinos too 



21 



22 



e+ u e Up , vT 



Vl 



e v e Vp. It 



23j: 7T+ — 

is easy to see that : Tf : Tf = 1/3 : 2/3 : 0. (Actually, v T are expected to be produced 
through the decay of Df charmed mesons generated also in pp and pry collisions. However, 



23|] and X* < HT 5 .) 



Df production is strongly suppressed 

The are very rapidly oscillating functions of energy. Taking into account the limited 
energy sensitivity of current and envisioned neutrino telescopes (AMANDA-II, for instance, 



had an energy resolution of 0.4 in the logarithm of the energy of the z/„-spawned muon 
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FIG. 1: (left) Eigenvalues 621 and 632 as functions of E , the energy at which the standard 
and NP energy-independent oscillation phases become comparable, i.e. </>y ~ Q° , according to 
Eq. (|18p . The redshift z = 1. The current upper bounds are plotted as horizontal lines. Notice 
that, due to these bounds, E NP cannot be lower than about 1 GeV. (right) Standard oscillation 
phase 021 and phase including the energy-independent contribution, </>2i , as functions of neutrino 
energy. The redshift z = 1. Note that the phases start to differ at E NF = 10 6 GeV, which 
corresponds to 621 = 6.1 x 1CP 29 GeV and 632 = 1.9 x 10~ 27 GeV. Below this energy, they are 
indistinguishable. 

we see that they are sensitive not to the instantaneous value of the ratios, (E ,z), but 
rather to the energy-averaged flavour ratios 



where E™ in = E Q - 5E a , E™ ax = E a + 5E D and AE Q = E" iax - Ef n = 25E D , with SE 



a small energy displacement. Without the NP effects, the high-energy neutrino flux from 
a distant astrophysical source is equally distributed among the three flavours, i.e. (T^) : 



In the presence of NP effects, however, the detected flavour ratios might be modified. 
Given that the relative dominance of the energy- independent NP phase over the standard 
phase 4>ij grows with energy, i.e. /<j>ij ~ E Q , we would expect that any modifications 
became more pronounced in the UHE range, PeV-EeV, or higher. As explained in Section 
HTl while the NP phase remains constant in energy, the standard phase decreases and, as a 
consequence, beyond a certain threshold (determined by the values of the the detected 
ratios would acquire a constant nonzero value, which might differ from the standard 




(17) 



<T£> : <T£) = 1/3 : 1/3 : 1/3. 



S 




E o [GeV] 



FIG. 2: Energy-averaged detected ratio T^, Eq. (fTT|) . as a function of neutrino energy E , 
for different values of the bij. Note that the ratio becomes constant only for unrealistically high 
energies: ~ 10 16 ' 5 GeV in the best case, when the bij are set to their upper bounds. For lower 
values of the bij, the energy at which the ratio becomes constant is higher. The neutrino flux from 
cosmic accelerators is predicted to span up to about 10 11 GeV; hence, the regime of constant 
due to an energy-independent contribution to the oscillation phase would not be observable. 



ratios 1/3 : 1/3 : 1/3, thus providing a distinct phenomenological signature of a possible 
energy-independent contribution to the oscillation phase. 

As a means of estimating the values of the for which the NP phase starts to be of 
importance, we can demand that Q° ~ <f>ij. From this requirement, we can calculate, for 
given values of the b^, the energy E NP above which the NP effects are expected to become 
increasingly more dominant in the oscillation. Doing this, we obtain 



_ g Am|.[eV 2 ]l-(l + 



-5/2 



E NP [GeV] = 6 x 10- 19 f^frw 1 - ^ ' ~' 3/2 . (18) 

bij [GeV] i_(i + ^)- 3 / 2 

The left panel of Fig. [1] shows a plot of 621 and 632 as functions of E . The current upper 
bounds are shown as horizontal lines. The lower the value of E NP , the earlier the NP effects 
would manifest. Notice that, due to the current bounds, E cannot be lower than about 
1 GeV. The plots have been generated for a fixed z — 1; for lower values of z we will have a 



higher value of by (30% if we take z = 0.03), while we will obtain a decrease in the values 
of bij for large z (a 20% decrease for z — 6). The right panel of Fig. [1] shows the standard 
and the energy- independent NP phases, 2 i and 2 i, respectively, as functions of neutrino 
energy, assuming that E NP = 10 6 GeV. Notice that the phases start to differ precisely at 
this energy. 

For concreteness, we will study the detected ratio of defined in Eq. (TIT)) , since our 
conclusions are independent of the chosen flavour. Fig. [2] shows the predicted ratio calculated 
for different values of &21 and 632 (we have assumed that 631 = 632 + &21) as a function of E . 

Fig. [2] shows that indeed becomes constant and different from 1/3 after a certain 
energy threshold. This occurs when the standard phase (j)y — > 0, so that, effectively, the 
oscillation phase is reduced to the energy-independent NP contribution, i.e. (py — > Q° , 
and the transition probabilities become constant. Note, however, that is constant 
only for E Q > 10 16 5 GeV in the most promising case, that is, when the by equal their 
current upper bounds. This is about five orders of magnitude higher than the energy of 
the most energetic neutrinos expected from cosmic accelerators. For smaller values of these 
parameters, the energy at which the ratio becomes constant is even higher. Using closer 
or more distant sources, effectively decreasing or increasing z, does not affect the energy 
threshold, but only modifies the constant value reached by T^. Therefore, we conclude 
that, given the current upper bounds on the by, an energy-independent NP contribution to 
neutrino oscillations would be visible in the high-energy astrophysical neutrino flux only if 
it modifies the oscillation amplitude (i.e. the mixing angles), as well as the phase. 

In light of this conclusion, within the formalism used in the present work, a comparative 
calculation, with and without NP effects, of high-energy astrophysical neutrinos detected at 
a second-generation neutrino telescope such as IceCube becomes unnecessary. 

IV. SUMMARY AND CONCLUSIONS 

We have considered the effect of a modified dispersion relation on the detected flavour 
ratios of high-energy neutrinos from cosmic accelerators. In the scenario of new physics 
that we have explored, the flavour oscillation phases are modified by the addition of energy- 
independent terms which depend on the parameters by. This contribution could correspond 
to a violation of CPT symmetry or to a nonuniversal neutrino coupling to a torsion field. 
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The current upper bounds on the 6^ are strict: 621 < 5.0 x 10 23 GeV and 632 < 1.6 x 10 21 
GeV. 

At sufficiently high energies, the oscillation phases are dominated by the energy- 
independent terms and the flavour ratios become constant and, possibly (depending on 
the values of the 6^ ) different from the average value of the ratios in the standard oscillation 
case, when new physics effects are absent. We have found, however, that even in the best 
case, when the bij are set to their upper bounds, the ratios are constant only for energies 
above 10 16 ' 5 GeV, about five orders of magnitude higher than the most energetic neutrinos 
that are expected from cosmic accelerators. Lower values of the will only result in higher 
energy thresholds for the ratios to become constant. 

Therefore, we conclude that, even though there could be, in principle, a clear signature of 
the presence of energy-independent contributions to the neutrino flavour oscillations, these 
are not detectable in the flavour ratios of high-energy neutrinos from cosmic accelerators if 
they affect solely the oscillation phases. 
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